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We investigated the flow profile of a polymer solution in a thinning capillary bridge. 
Fluorescent tracer particles with a diameter of 3/zm were used to visualize the flow. 
The cylindrical shape of the filament introduced strong optical abberations that could 
be corrected for, and we were able to characterize the flow in filaments with a thickness 
ranging from 150 to 30 fim. In the first regime when the filament was still sufficiently 
large, we used a PIV algorithm to deduce the flow field. At later stages when the 
number of particles in the observation plane decreased a PTV algorithm was used. 
The main two results of our measurements are as follows. First, the flow profile at the 
formation of the cylindrical filament is highly inhomogeneous and there is only flow 
in the outer parts of the filament. Second, we find that in most parts of the regime, 
where the temporal radius of the thinning filament can be fitted with an exponential 
law the flow indeed is purely extensional. 
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I. INTRODUCTION 



The Rayleigh-Plateau instability is the classical example for hydrodynamic instabilities in 
free surface flowi. It describes the disintegration of a liquid column into separate droplets via 
initial exponential growth of the sinusoidal deformation of the free surface. After sufficient 
viscous stretching, the system evolves into a nonlinear regime, and close to pinching, the 
system does not depend anymore on the primary configurations. The break-up process is 
described by self-similarity solutions. The addition of a small amount of polymers does 
not significantly alter the primary, linear instability because the solution has a very low 
linear elastic contribution and the polymers do not yet affect the flow at low shear rates 2 ; 
only, when the (elongational) flow rates become large enough the polymers intervene the 
flow. The thinning process is strongly slowed down, and instead of further necking of the 
capillary bridge, a cylindrical filament is formed that thins exponentially in time. The flow 
is supposed to be purely extensional. Figure CD, for example, shows the capillary breakup of 
a polymer solution. 




FIG. 1. Shadowgraph images from the capillary breakup of a 0.2 weight percent polyethylene oxide 
(PEO, Mw= 4 x 10 6 g/mol) in a 60/40 weight percent glycerol/water solution. Images are taken 
at t = -0.05, 0.25, 0.43 and 0.65 s (cf. Fig. 0}. The size of the images is 0.55 x 1.1 mm. The 
horizontal lines in the second image indicate the region shown in Fig. [5j The last image shows the 
final instability of the viscoelastic thread when many small droplets are formed. 



Extensional flow is of primary interest for the characterization of polymers because elastic 
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stresses can build up that are orders of magnitude larger compared to simple shear flow. 
The experimental determination of the elongational viscosity is still a challenging task and 
the characterization of this cylindrical thread is one of the rare reliable methods to extract 
at least some information on the elongational properties of polymer solutions. This exper- 
imental setup is also referred to as CaBER (capillary break-up extensional rheometer) and 
a commercial version is available as well. In spite of numerous experimental and theoretical 
investigations, there is still no complete understanding of the capillary break-up process^— . 
Most theories indeed predict an exponential thinning of the cylindrical thread^, but the 
determination of the correct time scales remains an open issued. Furthermore, the elonga- 
tional viscosity can be only deduced once a parallel filament has formed, but at this moment 
in time, the polymers have already accumulated a considerable amount of stress, and it 
is not clear how the transition from the primary thinning to the cylindrical filament shall 
be described. Numerical simulations indicate that there is a thin layer of large stress near 
the free surface of the thinning capillary bridge, while our experiments show the existence 
of a small layer with flow close to the free surface or no flow at all in the inner parts of 
the cylinder near the stagnation point where the flow is mostly elongational. Experimen- 
tal observations also show that the flow becomes very inhomogeneous at the final stage of 
the thinning process where the filament becomes unstable to the formation of secondary 
droplets and even some flow-induced phase separation has been reported^. However, our 
measurements show that in the intermediate regime, the flow indeed is purely elongational, 
and the only irregularity is the movement of the stagnation point. 



The experimental characterization of the flow field in a thinning capillary bridge implies 
certain difficulties due to the small sizes that even change in time and the optical aberrations 
due to the cylindrical form of the filament. In microscopic flows, fluorescent particles are 
often used to visualize the flow. Even if they are illuminated along the optical axes, the 
exiting light can be filtered before the camera, and only the fluorescence signal is observed 
that leads to a good contrast. The small sizes allow the use of microscopic objectives with 
a large numerical aperture. The desired plane of observation is determined by the field of 
depth, and there is no need to shape the incident light beam, e.g., in the form of a sheelr^ 1 ^. 

3 



II. THEORETICAL BACKGROUND 



A. Capillary breakup of a liquid bridge or droplet detachment 

The classical Rayleigh-Plateau instability describes the instability of a cylinder of infinite 
length, but the theoretical growth rate has been found to be in very good agreement with 
experimental data from a capillary thinning bridge 2 ' — . For an infinite thread, the instability 
evolves spatially in the form of a sinusoidal deformation, and for a short capillary bridge, 
the necking evolves in a symmetric manner. In the middle of the necking capillary bridge, a 
stagnation point is formed around which the flow is elongational. This elongational flow is 
very efficient in stretching the polymers, and when the flow rate times the polymer relaxation 
time, i.e. the Weissenberg number, is Wi = e\ > 0.5, the polymers intervene the flow. For 
such a flow, the Oldroyd-B model predicts a temporal evolution of the minimum filament 
diameter in the form: 

h(t) = /i e-* /{3Ao) . (1) 

ho is the diameter of the filament at the point when exponential thinning starts, and Ao is the 
polymeric relaxation time. The elongational viscosity can be deduced from the measurement 
of h(t) by the force balance of the capillary pressure and the polymeric stress. It follows 
that 

e is the elongation rate, and a the surface tension. A pure elongational flow in cylinder 
coordinates is defined by a velocity in axial direction and in radial direction in the form 

v z = ez (3) 

and 

v r = ~ir. (4) 

z and r are the axial and the radial position on the filament, respectively. Thus, the 
elongational rate can be determined by simply measuring the diameter h(t) of the filament 

m-^$. (5) 

4 



Together with Equ. [I] one finally gets an elongation rate which is constant in time: 

m = 3^ = e. (6) 

The exponential regime might last until the polymers are fully elongated. When the poly- 
mers inside the filament are completely stretched, the liquid might behave as a Newtonian 
fluid, but with a strongly increased viscosity. In this regime, the diameter of the filament 
is supposed to decrease linearly in time according to the self similarity solutions for viscous 
liquids 1 . 



III. EXPERIMENTAL METHODS 

A. Sample preparation 

The polymer in our study was polyethyleneoxide (PEO) with a molecular weight of My/ = 
4 x 10 6 g/mol (Sigma-Aldrich, Munich, Germany). PEO is a highly flexible polymer that 
has often been investigated in elongational measurements^^ 1 ^ 1 ^^— . The polymer was 
dissolved in a 60/40-weight percent glycerol-water solution at a concentration of c = 0.2 
weight percent, and the overlap concentration was c ov = 0.07 weight percent. Both the 
semidilute concentration^— and the viscosity of the solvent were chosen such that the 
characteristic time scales were slow enough to allow for a quantitative measurement of the 
flow profile. In general glycerol is hygroscopic, but this is only relevant at very concentrated 
glycerol solutions. The solution was prepared with the following protocol: the polymer 
was slowly poured in water first. After some minutes of swelling at rest, the polymer- 
water mixture was gently stirred for 24 hours. After that, the glycerol was added, and the 
solution was stirred gently again for 24 hours. The solution was measured within a few 
hours after preparation to minimize degradation of the polymers. This protocol lead to very 
reproducible results that were in agreement with data from the literature^ 1 ^. For the PIV 
and PTV measurements, a small amount of fluorescent particles (0.6 wt%) was added to 
the solution just before the measurement. The tracer particles had a diameter of 3/mi, and 
their absorption maximum was at a wavelength of 530 nm while the emission maximum 
was at 630 nm. The samples were characterized with a Thermo Haake MARS rheometer 
using cone plate geometries. The viscosity at a shear rate of 7 = 1 was rjx ~ 70mPas, and 
shear thinning was present up to a shear rate of 7 ~ 300 with 77 300 ~ 40mPas. Above this 
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shear rate, the fluid became unstable to a viscoelastic instability. The surface tension of the 
solvent was 7 ~ QhmN/m and its density p = 1150kg/m 3 . The addition of the polymer 
did not change the surface tension within the measurement resolution of the pendant drop 
method ±2mN/m. We should not expect that additional Marangoni stresses at the free 
surface of the viscoelastic filament affect the flow because the heat absorption of the laser 
was very low (< 1°C). We also tried to work with less concentrated or less viscous solutions 
but the dynamic was always to fast to allow for a quantitative PIV analysis. 

B. Experimental setup 

Our self-built Capillary Breakup Extensional Rheometer (CaBER) consisted of two cir- 
cular stainless steel plates with a diameter of 2 mm. A droplet of the polymer solution was 
put between the two plates. The upper plate was then moved away slowly with a speed 
of 0.05mm/s by a linear motor that was controlled by a PC. When the distance of the 
two plates overcame the critical length for Rayleigh-Plateau instability, the capillary bridge 
started to thin. This experimental protocol differs from the standard procedure where a 
steep step is performed to separate the plates and to introduce a kind of well-defined pre- 
stretch but similar to the slow retraction method presented in Ref.— . We find that our 
protocol is first more closely related to the more interesting case of a capillary bridge that 
thins due the Rayleigh-Plateau instability, and second, we find that the thinning of the 
filament is much more stable, e.g., there are less oscillations in the filament that disturb 
the PIV-measurements significantly. For the simple shadowgraph measurements (cf. fig- 
ure [1]), the thinning bridge was filmed with a high-speed camera that was equipped with a 
microscope objective with 5x magnification. 

For the flow visualization measurements, the fluorescent tracer particles in the solution 
could be illuminated either with a continuous Nd:YAG laser (A=532nm, P max ~250mW) or 
with a pulsed Nd:YLF laser (A=527nm, P m ax ~16W). We chose a transmission geometry, 
and the part of the laser light that was not diffracted by the filament passed through the 
sample before entering the 20 x magnification microscope objective with an numerical aper- 
ture of NA = 0.4 (FigfJJ). A filter cube that was located between the objective and the tube 
lens (see Fig. [2]) contained a 605/55nm band-pass filter that protected the camera chip from 
the laser radiation, but let pass the fluorescent light. 
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FIG. 2. Sketch of the experimental setup. 
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FIG. 3. Inset: Cylindrical lens on millimeter paper. Squares: Measured radial positions of the 
vertical lines recorded with and without the cylindrical lens on the millimeter paper. The line is 
given by Equ. [7J 



Due to the cylindrical shape of the filament, the image of the fluorescent tracer particles 
is distorted. We only measured in the middle plane of the filament, and the effect can be 
visualized with a half-cylinder lens. Figure |3] shows a half cylinder lens lying on a sheet 
of millimeter paper. Plotting the radial positions of the lines recorded via a cylindrical 
lens versus those without a lens, one gets a straight line, at least if one is not too close to 
the radius r. This is in good agreement with a theoretical consideration based on simple 
geometrical optics. A point source at a distance x from the center of the lens perpendicular 
to the optical axes seems to be shifted by a length Ax that is given by 
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Ax = rcos yarcsin ^ an (arcsin (n^—j arcsin i^j^j ■ (?) 

n 2 is the refractive index of the cylinder lens. Equ. [7]predicts a close to linear dependency 
of Ax(x) if x is smaller than x 3> r/n 2 . This relation allows for a simple correction of the 
geometrical aberrations. However, we were still limited to the regime in which the parallel 
filament already had been formed. Before the formation of the parallel filament, there is an 
additional curvature in the direction of the longitudinal axes that is more difficult to correct 
for. 

In the course of thinning of the filament, the number of observable particles was not 
constant, and we had to start with a comparably high concentration of particles of 0.6 weight 
percent. A PIV analysis was only possible at the beginning of the thinning process. During 
the breakup, the number of observable particles decreased, and a PTV analysis was more 
reasonable from the moment when single particles were distinguishable. Devasenathipathy 
et al.— used /iPIV as well as PTV to characterize the flow profile at the intersection of a 
cross-channel. Their measurements showed that the two techniques provide the same results. 
This was also confirmed by simulations and is in accordance with our findings. 

The images were captured at a rate of up to 150 Hz and afterwards analyzed by the 
MatPIV algorithm^. In our measurement, a long but thin filament was investigated, and 
the interrogation window width was set to 16 pixels in the radial direction and 64 pixels in the 
axial direction. The overlap of the windows was 50% for both directions. The depth of the 
interrogation plane 5 was determined by the numerical aperture of our microscope objective, 
and for the 3/xm particles we get 5 = 2A/j,mr&. The interface between the areas with and 
without particles moves in the radial direction because the filament thins in time. This 
movement might lead to spurious velocity vectors in the analysis such that non-vanishing 
velocity components outside the filament can appear that have to be sorted out manually 
or by appropriate filtering. For the PTV analysis the centroid algorithm from Lab- View 
(National Instruments, Austin, Texas) was used. 

The smallest capillary threads that were evaluated had a diameter of h ~ 30 fim and the 
concentration of the microbeads were low enough to not expect any alteration of the flow due 
their presence, similar to what is found in standard microfluidic experiments^. Furthermore, 
we did not observe any concentration enhancement of the microbeads near the surface of 
the viscoelastic thread and we assume that the surface tension is not significantly affected 



by the microbeads. 

IV. MEASUREMENTS AND RESULTS 
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FIG. 4. Filament diameter versus time. The dashed line is taken from the simple shadowgraph 
images. It shows the beginning of the capillary instability before the parallel filament has formed. 
The squares are taken from the background signal of the PIV images (Fig. [5]). The straight line 
depicts an exponential fit to the experimental data in the range 0.1 < t < 0.5. The characteristic 
time scale Ac = 227ms is calculated according to Equ. HJ an d the elongation rate is e = 8.81s _1 . 
The black arrows below the curve (ticketed with (1), (2) and (3)) mark the moments at which the 
pictures in Fig. [5] are captured. The two left arrows above the curve (PIV) label the transitions 
in the PIV analysis, which are shown in Fig. [6] and in Fig. [8l The two right arrows above the 
curve (PTV) locate the temporal borders in between which the PTV analysis shown in Fig. [121 is 
performed. 

A. PIV analysis 

An overview over the thinning process is given in Fig. HJ It shows the filament diameter 
as a function of time and indicates at which time the data are taken. Figure |5] shows three 
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PIV-images that are captured at different moments during the capillary breakup of the 
polymer solution. Obviously, the number of particles strongly decreased in the course of 
time. Still, the filament was shining fluorescently everywhere because there was not only 
fluorescence light from the particles in the focal point but also from particles that are out 
of the focal plane. This effect allowed the determination of the radius h(t) of the filament 
and verification of whether the plane of measurement was in the middle of the filament. For 
the first two pictures, only PIV analysis was possible, but for the third picture, both a PIV 
analysis and a PTV analysis could be performed. 

The axial position of the stagnation point of the elongation flow in the first image of Fig. 
[5] is at a hight of about 150/xm. Its exact position is difficult to determine, but apparently 
it was moving upwards and it has vanished from the image section in the third picture. 
According to Equ. EJ the velocity v z at a fixed axial position should be constant for all 
times because the elongation rate is a constant. In Fig. El apparently v z at a given position 
is not constant in time. The flow has not yet completely formed in the first images, and the 
neutral (stagnation) point moves in time. For a more quantitative analysis, a fixed axial 
position z at z = 600fim was chosen, and the axial velocity v z as a function of the radial 
position r was determined. The result is shown in Fig. [6] for three profiles at the beginning 
of the breakup. Positive values mean that the flow goes downwards, whereas negative values 
represent an upward flow. The three profiles show that there is only a flow near the edge of 
the filament but no flow in the middle. For t=10 ms and t=20 ms, the liquid flows upwards, 
whereas, at t=30ms, it flows downwards. The axial velocities are still small, but the flow has 
changed its direction between 20 and 30ms after the beginning of the measurement. This 
is because the neutral (stagnation) point has moved upwards as well, and the stationary 
measurement level is first below and then above the stagnation point. 

An inhomogeneity in the flow profile along the r-axes could be explained by at least two 
different mechanisms. First, the amount of stretch at the boundary is largest due to the 
larger deformations compared to the center of the filament. This should lead to a stress 
layer near the boundary. This seems to be indeed the case in some numerical simulations^. 
Another possible explanation might be that the flow close to the center and the stagnation 
point is closer to a pure elongational flow that is more efficient in stretching the polymers 
and causes a larger elastic stress that suppresses any flow. This could be closely related to 
what is known from stagnation point flow and is referred to as the " birefringent strand"—. 
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FIG. 5. PIV images at three different moments during the capillary breakup (see FigH]). At the 
lower part of the filament a fully developed elongational flow is present. Closer to the stagnation 
point in the upper part of the filament there is stronger axial flow at the boundaries then at the 
center of the filament. The axial position of the stagnation point moves upward in time. Please 
remark that the arrows are rescaled from image to image. 



It refers to a region of high birefringence (and thus high elastic stress) along the center line 
close to a stagnation point. 

With increasing time, this profile changes qualitatively (Fig. [7j). The maximum of the 
axial velocity increases from 0.1 mm/s up to about 3.5mm/s and then remains constant. 
Moreover, the spatial inhomogeneous profile with some flow at the edge and no flow in the 
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FIG. 6. PIV profiles of the axial velocity v z at the beginning of the measurement between t=10ms 
and 30ms. In this time interval, the filament diameter changes from h = 168 to 143/J.m. 




FIG. 7. PIV profiles of the axial velocity v z between 10ms and 250ms with a time-lag of 30ms. In 
this time interval the thickness of the filament changes from h = 168/im to 48/^m. 

middle of the filament that has formed at the beginning (cf. Fig. [6]) is conserved only 
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for some time. First, the profile grows until the velocity at the edge has reached a value 
of about 2.5 mm/s. Then, the expected flow with a nearly constant velocity for all radial 
positions at a fixed axial position is formed such that the flow "grows" continuously into 
the center of the filament (Fig. [HD- This growth occurs in a manner that is not fully 
symmetric, but rather growth from one side is preferred. The change from a zero velocity 
to the finite maximum value at a position r next to the transition from no-flow to flow 
occurs in a short time interval, i.e., within the 6.5 ms measurement resolution. This could 
compared to the time scale of viscous diffusion, which is even much shorter considering 
the high extensional viscosity present here r = h 2 /(r) e /p) ~ lfis. The complete transition 
from the inhomogeneous profile lasts approximately 20 ms after the boundary velocity has 
reached its maximum at t — 170ms after the beginning of the measurement. 



2.5 

2.0 

^ 1.5 
w 

E 

J5, i.o 

N 

> 

0.5 
0.0 

-40 -30 -20 -10 10 20 30 40 50 60 

radial position (pm) 

FIG. 8. PIV profiles of the axial velocity v z between 150 ms and 170 ms that show the transition 
from the pure edge flow to the plug flow-like profile. 

If we evaluate the velocities along the filament (Fig. |9]) we see that there is a region 
without flow and only at some distance from the stagnation point there is flow. The velocity 
profile in Fig. [9] already suggests that the flow is purely elongational in this regime. The 
velocity v z along the vertical axis can be fitted with the linear law Equ. [3j The slope of 
this linear fit directly gives e. Doing so for the central layer of the filament (parallel to the 
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FIG. 9. Velocity along the centerline for different times. The elongational flow implies a linear 
dependency of the axial velocity on the axial position. The slope of the linear fit give the elonga- 
tional rate e (see Fig. I1U|) . Close to the stagnation point, there is a region without any flow and 
this region moves upward in time (see Fig. [TTJ) 

z-axis) for different times during the thinning process, one gets the results shown in Fig. [TQl 
The straight line marks the value of e that is obtained from the relaxation time Ac (cf. Fig. 
H|) that can be deduced from the measurement of the filament diameter h(t) by use of Equ. 
[6j There is a sharp transition at t = 160ms. We do not have an explanation for this sharp 
transition, but this is roughly the moment when there is flow in most part of the filament. 
The data at times t < 160ms are taken when the flow profile all over the filament is still 
inhomogeneous; in principle, one can not yet define a constant elongational rate over the 
whole filament. For a typical CaBER application, the exponential fit of the diameter h(t) 
is the most important quantity, and apparently one gets reasonable results if one chooses 
the range over which the diameter h(t) thins exponentially by visual inspection (here we 
choose 100ms < t < 500ms). If, for example, one chooses the starting point of the time 
interval by considering the analysis above and shifts the end for the same amount such that 
170ms < t < 430ms, the characteristic time scale changes only by 0.5%. 

In Fig. [11] we plot the position of the neutral point and the position of the transition 
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FIG. 10. Elongation rate at the center of the filament. Squares are calculated by linearly fitting 
the experimental PIV data of v z (Fig. [9|) with respect to z according to equation [3] . The straight 
line marks the value of i that is obtained from the exponential fit of the diameter (Fig. UJ. 
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FIG. 11. Position of the stagnation point (squares) and position of the region along the center line 
with the transition from the elongational flow to zero flow (triangles, see Figl9]). 
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regime from no flow to flow along the centerline. The position of the neutral point is not 
defined a priori and in principle it can move freely. Though this implies acceleration of fluid 
elements, the stagnation point indeed moves upward to the upper half droplet. With this 
movement, the transition regime of no flow to flow moves upwards too and therefore the 
fully developed elongational flow evolves only gradually at a given position. 

B. PTV analysis 

From t = 160ms the flow should be purely elongational and the velocity v has both an 
axial and a radial component. We tried to deduce the radial component v r from our PIV 
data, but v r <^v z , and we were not able to resolve v r in a quantitative manner. In Fig. 
the arrows are nearly parallel to the z-axis, and finally, we had to use a PTV analysis to 
quantify both velocity components. The PTV analysis could be used only at longer times 
when the number of tracer particles became small enough. Furthermore, we wanted to verify 
if the streamlines are in accordance with the assumption of a pure elongational flow. The 
result is shown in Fig. [12] and [T3j Resolving Equ. [3] and H] one finds by integration: 

where zq and tq are the initial values in the axial and radial direction, respectively. 
However, in our coordinate system the trajectories do not start at (r — > ±oo,z — > 0) and 
end at (r — > 0,z — > ±oo), and the experimental trajectories are given by 



Zq again is the initial value in the axial direction (for r — > ±oo), and is the final value 
in radial direction (for z — > ±oo) . The latter value is zero for particles on the left side of the 
filament and positive for those on the right side. The factor a is a free fit parameter that is 
always positive because z is the minimum value in axial direction, and (r — r^) 2 > 0. 

Both trajectories of a particle on the left- and right-hand sides of the thinning filament 
are shown in Fig. [131 an d a fit according to Equ. [9] matches the experimental data very 
well. 



a 



(9) 
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FIG. 12. Trajectories of different single particles determined by the PTV analysis. The open 
symbols are the trajectories that are shown in Fig. PT3l 
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FIG. 13. Trajectories of a particle on the a) left-hand side and on the b) right-hand side of the 
filament determined by the PTV analysis. The solid line is a fit according to Equ. [9j 

V. CONCLUSION 

In conclusion, we presented PIV and PTV data of a thinning capillary bridge of a polymer 
solution. The optical aberrations that are introduced by the cylindrical filament could 



17 



be corrected for, and the thinning process could be characterized for a filament thickness 
reaching from 30/im < h(t) < 150/im. The measurements at the beginning of the breakup 
that are analyzed via PIV show that the axial velocity v z is first strongly varying. Initially, 
the liquid has a finite velocity only near the edge of the filament, whereas there is no flow in 
the middle. A possible explanation could be that, before the formation of the filament, there 
is a stagnation point in the center of the capillary bridge and the stretching and build-up 
of stresses are most efficient close to the center. Both the slight asymmetric flow profile 
apparent in Fig. [7] as well as the non stationarity of the stagnation point might be caused 
by gravitational sagging. The Bond number Bo = pgh 2 jo ~ 10~ 3 is very small, but former 
studies on the final disintegration process of the viscoelastic thread showed that gradational 
effects still play a role e.g. for the formation of the differently shaped transition regions from 
the filament to the upper and lower half droplet at the ends of the filament 31 . 

Later on, the inhomogeneous profile evolves into a homogenous one, and the elongation 
rates (at a fixed axial position) are constant over time for each radial position. The radial 
velocity v r was found to be much smaller than v z , and it was not possible to characterize v r 
quantitatively with PIV. Therefore, we used PTV to extract the flow lines, and we were able 
to show that the flow, once fully developed, indeed is a pure elongational one. This happens 
approximately at the moment when an exponential fit of the filament diameter h(t) becomes 
reasonable. This quantity is used in simple CaBER (capillary breakup extensional rheology) 
experiments to deduce the elongational rate and more importantly the elongational viscosity. 
For such a type of material characterization, we could show that the flow indeed is purely 
elongational, at least for our polymeric system. 
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